Structure Factor in the Presence of Shear - an RPA Calculation 
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o : 

, We consider the structure factor of a system of colloidal particles immersed in a host liquid. Each 

fSJ . particle is assumed to be affected by forces due to other particles, a drag force proportional to the 

» I ' velocity of the particle relative to the local velocity of the fluid and a fluctuating random noise. The 

CL(, effect of the particles on the velocity field of the host liquid is neglected. The problem is treated 

■ within the RPA approximation which is a widely used tool in many body theory. The validity of 
the RPA result is discussed. 
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, PACS Numbers: 83.50.A, 82.70.Dd, 64.60.My, 61.20.Ne 



C/2 



We consider a system of identical particles immersed in a host liquid subject to an external linear shear flow. Such 
systems have beed studied experimentally and numerically for the last two decades - 0- The results of those 
O ' studies showed not only distortion of the structure factor due to the external shear, but also more dramatic effects like 
' shear induced ordering and jamming. The theoretical explanations were first given in a paper by S. Hess H], followed 
I • by a number of other papers that tackled the problem from different angles - [ p^ . It seems, however, that in spite 
[ of substantial advances made over the years, a simple and full theoretical description is still lacking. In the present 
paper we will obtain the shear dependent structure factor using an RPA approximation employed on an extension of 
the exact collective coordinates Fokker-Planck equation derived recently by Edwards and Schwartz jl9j. The whole 
d , method is based on borrowing approximations that proved simple and effective in the study of quantum mechanical 
^ ' many body systems. Such methods have been applied in the past |20), and much more simply recently to 
I ' obtain the Percus-Yevick (PY) equation for a classical hard sphere system. We hope to produce the analog of 
■ O the hard sphere PY equation for the case of shear flow in the very near future. 

^ I Our starting point is an assumption concerning the forces acting on each particle. We assume that each particle 
is affected due to forces applied by the other particles, by a drag force proportional to its velocity relative to the 
local velocity of the liquid and by a fluctuating noise which is assumed to be uncorrelated for different particles. It is 
further assumed that the velocity field of the liquid is given externally and is unaffected by the motion of the particles. 
This assumption will be checked in a future publictaion. We can say at present, that the main effect we expect from 
the inclusion of the modification of the velocity field by the system of particles is a modification of the interaction 
between particles and a modification of the noise. In any case, our present model is described by the set of Langevin 
' equations 

^ ■ l[r^ - V{n)] = -V, )_J win - rj) + m{t), (1) 

o ■ 

d ' where n is the radius vector of the «'th particle, V{r) is the divergenceless velocity field of the liquid, 7 is the friction 
^ , constant, N is the number of particles, w is the two body potential between particles and f]i{t) is a random force 
acting on the i'th particle obeying 

g: (^fw) = o (2) 

O ' and 
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5^ ■ where k and I denote cartesian components of a vector. The Langevin equation above leads in a standard way to a 
^ Fokker-Planck equation for the probability distribution of particle coordinates, P 



dP _ d 
dt 7 ^ dr^k 



dP dU 
kT-^ + -^P~jVk{n)P 



^ LP, (4) 
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is the total potential energy, and kT 
functionals of the density 



Since the particles are identical, the only relevant physical observables are 



N 



(5) 



To be specific, we assume that our A'^ particles are enclosed in a cubic box of volume V with periodic boundary 
conditions. The natural collective coordinates are the Fourier transforms of the density 



p^{{n}) = N-^ I p{f)e-'^-'d^r for with = 

' V 3 



(6) 



The Fokker-Planck equation for the probability distribution of particle coordinates can be transformed into a Fokker- 
Planck equation for the probability distribution of collective coordinates by following the procedure outlined by 
Edwards and Schwartz We note first, that the probability distribution of the collective coordinates, P{{p^}), is 
given by 



p{{p^}) ^ I p({ri}) n ( - E I n 



(Pri. 
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Therefore if 0({r^}) is an observable that can be written as 0({/5g({?^})}), then its thermodynamic average can be 
calculated by 



(0) 



(8) 



(A similar transformation was used previously for the case of equilibrium statistical physics ) . We now obtain the 
Fokker-Planck equation for P{{pq}) by multiplying eq. (^) by 



N 



-iq.rj 



and integrating over the particle coordinates. In ref. jl^] that considered a system without shear, an explicit equation 
is obtained for P{{p^}). Therefore, here we need only to introduce the term corresponding to the shear. Since the 
Fokker-Planck equation is linear, the form we obtain is 



LoP + iN-^y^Qk^ 



dt 



I e-^«^-^Vfc(rl)P({r;„}) J] 5 [pp-N-^- ^ e^^^"^ | \{d 

•' P#0 \ 3 = 1 I 1=1 



n 
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where Lq is the linear operator obtained by Edwards and Schwartz ||T^. We do not dwell here on bringing the term 
due to the shear for a general divergenceless velocity field to the form ALP. We consider instead the special case of 
linear shear 



V{f) = Cxz. 

In that case, the last term on the R.H.S of eq. becomes 

^ dp,-dP 

^ dqi dpff 

q 

The full form of eq. (||), including the specific Lq is thus given by 



(10) 



(11) 



dP kT 



dt 
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dp^ dP 
dqi dpil 
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d^P , P ,^ dP 



(12) 
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where v = ?S and ^^{q) is the Fourier transform of the inter-particle potential defined by 



(13) 

We will look for a stationary state of the system, hence the L.H.S vanishes. In order to keep the discussion simple, we 
will employ the RPA approximation on the R.H.S of eq. ( p^ ) by discarding all the third order terms except for those 
where one of the p's is po — VN. (Note that derivatives with respect to po do not appear). Thus, we're left with the 
following equation for the stationary probability distribution: 



E 



A-\q)P + A-\q)p,^Z- + 



where A ^{q) = 1 



kT 



W 



"^dpif ' dpg'dp^q- 
{q). We try a solution of eq. (|lj) of the form 
P({p,^}) = e-^S?^o«^'(?>f''- 



' dqi dpff 



(14) 



(15) 



(We denote the pre factor of PqP-q in the above Gaussian form by S^^{q) since if such a form does realy solve eq. (|TJ) 
then that pre factor will be, as follows from eq. (||), the inverse of the structure factor we are looking for.) The 
combination of equations (^ and ([l^ ) yields the following equation for the structure factor: 



- S-\q) - A-\q)S-\q)pq^p^q-+ S-'{q)Pq^p^^] + „q,S-\q)^^p^^'j 



0. 



(16) 



The last term in eq. (16) can be rewritten as 



« 9 



Integrating eq. (|17l ) by parts and combining the result with the rest of the terms in eq. ( |16D gives 

dSiq) 



Y^lq' [A-\q-) - S-\q-)] 



-q'A-\q)S-\q) + q'S'^q) + -q^S-^iq)- „ 

2 dqi 



PqP-q 



(17) 
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This equation for the static structure factor should hold for every collective coordinates configuration, {p^}, hence 
the following conditions must be obeyed: 



Y^q^ [A-\q)-S-\q)\ =0 



(19) 



and 



5(g) vq^dSjq) 
A{q) ^ 2 g2 dqi 



= 0. 



(20) 



Although it seems that the fact that there is a redundant condition here may cause trouble, it can be easily verified 
that any solution of the second condition also satisfies the first condition. Viewing the second condition as an ordinary 
differential equation in qi leads to a solution which depends on an initial condition. This condition is determined by 
applying the constraint that the structure factor must be finite for all the s. The structure factor turns out to be 



S{q) = — 



(21) 



where (f = (9^,(72,93) and q" = (q" , 92, ga)- This is the RPA dependence of the structure factor on the shear rate. 
A number of points have to be addressed now. The structure factor is given here in terms of the Fourier transform 
of the inter-particle potential, w{q}. In the case when w{q) does exist, one may replace A~^{q) in eq. ( ^l| ) by the 
simplest RPA approximation for the structure factor in the absence of shear 
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S^\q) = \ + ^w{q). (22) 

Indeed, in the limit of zero shear rate, our calculation does yield the well known RPA result (|2^). However, The 
experience with non sheared liquids suggests that even when w{q) exists, in most cases such RPA approximations 



produce bad results. A useful alternative is to turn eq. (22) around and to calculate an effective potential out of the 



externally given structure factor. This procedure actually replaces w(r) by the Ornstein-Zernike direct correlation 
function. Rewriting eq. (|l]) using an externally given S^^ici) recovers exactly the expression given by Ronis for 
the structure factor in the presence of shear. However, a further study of the Ronis expression reveals an interesting 
discrepancy when hard spheres interactions are considered. The structure factor in the absence of shear, So{q), must 
have the property that the associated pair distribution function vanishes within the hard spere diameter. Further 
more, S{cj), the structure factor in the presence of shear, must also have this property. Unfortunately, it seems that 
this property does not follow from the Ronis expression when inesrting the external expression for Sq (q) given by 
Wertheim ^ - Thiclc [||| into the Ronis expression. We expect to deal with this problem in the very near future. 
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